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ATTITUDE DETERMINATION OF THE SPIN- STABILIZED 
PROJECT SCANNER SPACECRAFT 

By T. M. Walsh, Jean C. Keating, 
and Dwayne E. Hinton 
Langley Research Center 

SUMMARY 

The Scanner experiment required continuous knowledge of the spacecraft celestial 
attitude with a one-sigma accuracy of 0.016°. The basic measurements used to obtain 
attitude information were sighting times of stars viewed by the star mapper on the Project 
Scanner vehicle. Each star sighting produced coded-pulse sequences the amplitudes of 
which were related to the visual magnitude of the star. The times of occurrences of these 
pulse sequences were related to the angular separation between a pair of stars. These 
angular-separation data were used to identify the sighted stars. 

The right ascension and declination of the identified stars and their transit times 
were used to identify the values of a set of parameters in an assumed model of the vehi- 
cle dynamical equations of motion. These vehicle motion parameters were then used to 
determine the vehicle celestial attitude as a function of time. 

A brief discussion of the star mapper is presented and the technique used to detect 
transit times of stars from the coded signals of the star mapper is described. A complete 
description of the star identification and parameter identification procedures is presented. 
A preflight analysis, used to establish predicted accuracies of attitude determination, is 
described. The estimated accuracies of vehicle attitude and star- mapper optical- axis 
pointing direction for Project Scanner flights 1 and 2 are discussed. The estimated one- 
sigma error in pointing direction of the star-mapper optical axis for these two flights was 
0.006°; this value compared favorably with the experimental requirement of 0.016°. 

INTRODUCTION 

Langley Research Center has been engaged in a series of experiments designed to 
obtain flight measurements of the infrared radiance signature of the earth's horizon. 

(See ref. 1.) The intent of these experiments is to provide necessary basic information 
to design and evaluate infrared horizon sensors which utilize the earth’s limb or horizon 
for spacecraft attitude-measurement purposes. Project Scanner was designed to obtain 
high resolution measurements of the earth's horizon radiance profile in two infrared 


spectral intervals and to position these measurements accurately with respect to the solid 
earth. The Scanner experiment was flown on a suborbital ballistic rocket vehicle fired 
from Wallops Island, Virginia. The first flight was made on August 16, 1966, and the 
second, on December 10, 1966. 

An operational schematic of the Scanner spacecraft is shown in figure 1. The space- 
craft was despun to approximately 3/4 revolution per second when horizon data-gathering 
altitude was reached, and was erected to within 2° of the local vertical by a cold gas reac- 
tion jet control system similar to that described in reference 2. Two radiometers, each 
having a basic resolution of 0.025°, were mounted back to back in the spacecraft. Radiom- 
eter 1 was used to sense radiance in the 14- to 16-micron or carbon dioxide region and 
radiometer 2 was used to sense radiance in the 20- to 35-micron or water vapor region. 
The fields of view of the radiometers were scanned through the horizon at 10° per second 
by motor-driven mirrors. Precise attitude and submissile point information was needed 
to position the infrared radiance measurements at earth geographical locations. The geo- 
graphical position of the spacecraft submissile point was determined with the aid of stan- 
dard radar tracking techniques. Celestial attitude information was obtained by means of 
the star mapper in the Scanner vehicle (ref. 3) which was used to detect star transits as 
the star-mapper field of view was scanned through the celestial sphere by the spacecraft 
spin. All spacecraft data were transmitted to earth on a real-time basis for data pro- 
cessing and later combined with ground tracking data for reconstruction of the measured 
earth infrared horizon profiles. This report is concerned primarily with a description of 
the data obtained from the star mapper and the data processing necessary to determine 
the spacecraft celestial attitude. Related discussions are presented to show the develop- 
ment of the data- reduction accuracy criteria and the results of application of the data pro- 
cessing to flight data. 


SYMBOLS 


A,B,C transformation matrices 

a estimated azimuthal separation of a pair of observed stars, degrees 

di,d 2 earth radii, kilometers 

estimated elevation angle of kth observed star in star- mapper field of view, 
degrees 

F a matrix of first partial derivatives of JDEHR with respect to each undeter- 

mined element of X 


2 


an e i em ent of matrix F 

G 1> G 2 any two stars of celestial sphere (fig. 11) 

g c angular separation of a pair of stars based on cataloged star positions, 

degrees 


g Q estimated angular separation of a pair of stars based on star- mapper data, 

degrees 


Ag 

h 

h t 

Il,l2,l3 

A A 

i,j,k 

J,D,E,H 

M 

m v 

<4 

R 

r 

T l 

T 

t 


allowable difference in angular separation of cataloged stars and estimated 
angular separation of observed stars, degrees 

spacecraft altitude, kilometers 

tangent height, kilometers 

spacecraft moments of inertia, slug-ft 2 

unit vectors with subscripts denoting specific frames (fig. 6) 

transformation matrices (table HI) 

tolerance limit on convergence of correction vector AX 
apparent visual magnitude of a star 

jth estimate of quadratic form of error in k star sightings based on vec- 
tor r 

unit vector to a star in celestial reference frame 
unit vector to a star in a slit plane reference frame 
an element of unit vector r 

time duration of a star transit across narrowest reticle slit, seconds 
time, seconds 
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tSk^Vk 

4 

At k 

V 

X 

X j 

AX 

AX j 


x m 


Ax m 
AY 
Ay 
a, 6 
0 


01,02 


0s>0v 


observed kth star transit time, seconds 

observed star transit time of slanted and vertical slit planes, seconds 
true kth star transit time, seconds 
a specific instant of time, seconds 
kth time residual, seconds 

matrix multiplier of modified vector of undetermined parameters, AY 

vector of undetermined parameters 

jth estimate of vector of undetermined parameters 

correction vector of undetermined parameters 

jth estimate of AX 

an element of vector of m undetermined parameters 

an element of correction vector AX 

modified correction vector of undetermined parameters 

element of modified correction vector AY 

right ascension and declination, respectively, degrees 

rotation of a slit plane with respect to a plane perpendicular to an arbitrary 
azimuthal plane, degrees 

rotation of slanted and vertical slit planes, respectively, with respect to a 
plane perpendicular to test table azimuthal plane, degrees 

rotation of slanted and vertical slit planes, respectively, with respect to a 
plane perpendicular to star-mapper azimuthal plane, degrees 
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y azimuthal displacement of any slit plane from optical axis, degrees 

7S>>V azimuthal displacement of slanted and vertical slit planes from optical axis, 

degrees 

e l’ e 2 misalinement angles of vehicle principal spin axis relative to star- mapper 

reference frame, degrees 

e tk simulated kth star transit time error, seconds 

time hi story of pointing-direction error of vehicle principal spin axis in 
right ascension and declination, respectively, degrees 

e a (t)8,e 6 (t)8 time history of pointing-direction error of optical axis in right ascension 
and declination, respectively, degrees 

V star elevation angle in any slit plane, degrees 

Vs^y star elevation angle in slanted and vertical slit planes, respectively, degrees 

6,(f),il/ Euler angles defining orientation of vehicle principal axis relative to an 

inertial reference frame, degrees 

0,$ Euler angles defining the position of vehicle angular momentum vector in 

celestial reference frame, degrees 

M angular width of narrowest reticle slit, degrees 

°a standard deviation of pointing-direction error of any selected axis, degrees 

CT a 6 ? °a(t is ) 8 simulation results of standard deviation of pointing— direction error of 
vehicle principal axis and optical axis at specific instants of time, 
degrees 

a a^^)max maximum standard deviation of pointing-direction error of optical axis 

at specific instant of time for simulated data, degrees 

°t standard deviation of time errors used with simulated data, seconds 
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°tr 


standard deviation of time residuals for simulated or flight data, seconds 


standard deviation of pointing-direction error of vehicle principal axis in 
right ascension for a specific instant of time of simulated data, degrees 


standard deviation of pointing-direction error of optical axis in right ascen- 
sion for a specific instant of time of simulated data, degrees 


<p 0 ,\p 0 initial values of <p and ip, degrees 

oj total spin rate of vehicle principal spin axis, degrees/second 

w t estimate of u> from raw star-mapper signals, degrees/second 


Subscript: 


k integer 


Superscripts: 

T matrix transpose 


-1 matrix inverse 

Derivatives with respect to time are denoted by dots over the variable. An under- 
line denotes a vector, a ~ denotes an estimated value and ~ denotes a unit vector. 

PROBLEM DESCRIPTION 

The experimental requirement of locating the earth’s horizon radiance profiles rela- 
tive to the solid earth can be best expressed in terms of locating a radiometer line of sight 
relative to an earth radius. This relationship is shown in figure 2 where the radiometer 
line of sight is shown displaced from the solid earth by a distance ht defined as the tan- 
gent height. The spacecraft is indicated at an altitude h, and the earth's radius is indi- 
cated by the radii di and d 2 . The Scanner experiment required that the measured 
radiance profiles be located relative to the solid earth with a one-sigma accuracy of 
1.5 kilometers for a spacecraft altitude of 1100 kilometers. This accuracy requirement 
may be expressed in terms of an angular accuracy in positioning of the radiometer line of 
sight as 0.022°. The significant error sources which affect this accuracy are: errors 
in radar tracking of the spacecraft altitude and submissile point, uncertainties in the 
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radiometer scan mirror position measurements, uncertainties in the knowledge of the 
alinement between the star mapper and radiometer, and inaccuracies in the celestial atti- 
tude determination from using star- mapper information. The contributions of all error 
sources excluding celestial attitude errors were estimated to be 0.015°. The allowable 
tangent height error due to errors in celestial-attitude determination was determined by 
assuming that all tangent height errors were independent and that their root-sum-square 
value was equal to the total allowable error expressed in degrees. By using this assump- 
tion, the allowable one -sigma error in spacecraft celestial attitude determination by using 
star-mapper information was computed to be 0.016°. Past experience with the difficulties 
in achieving the theoretical design limits of attitude -measurement systems by using com- 
plex optical devices led to the more conservative design goal of an allowable one-sigma 
error of 0.006° in spacecraft celestial attitude determination. 

The measurements used to obtain attitude information were provided by a passively 
scanning star telescope called a star mapper, which is described in detail in reference 3. 
The star mapper is shown mounted in the spacecraft in figure 1. The instrument used a 
refractive optical system based on a modified Petzval design with a 15- inch (38.1 cm) 
focal length, a 5-inch (12.7 cm) clear aperture, and a field of view of 6° by 6°. The 
instrument- mounting positioned the star-mapper optical axis normal to the expected prin- 
cipal spin axis of the spacecraft. A coded reticle, shown in figure 3, was centered in the 
focal plane of the instrument, and a photo detector was placed behind the reticle. The 
reticle was principally opaque and contained a pair of identically coded sequences of 
transparent slits, one vertical group and one slanted group. The narrowest slit had a 
resolution of 0.015°. The coding feature of the reticle was used to reduce the system 
response to noise and to the multitude of stars dimmer than a threshold visual magnitude; 
this feature was selected to insure a minimum of three detectable and identifiable stars 
per vehicle spin period. (See ref. 4.) Spin motion of the spacecraft caused star images 
to travel across the reticle and produced two coded sequences of radiant energy which 
were sensed by the photo detector and transmitted to a ground receiving station. The 
ground-received coded star signals were tape recorded simultaneously with range time 
signals having a 10-microsecond resolution. 

The procedures involved in determining celestial attitude of the Scanner vehicle 
from this raw data tape are illustrated in figure 4. Processing of the recorded star sig- 
nals consisted of four major operations: time detection of apparent star sightings; manual 
pairing of star sighting times to reduce the number of false star sightings; identification 
of detected stars; and identification of the parameters of the dynamical equations of motion 
describing the spacecraft attitude. A brief description of the first two operations, which 
are discussed in reference 3, is presented herein for clarity, and the latter two operations 
are discussed in detail in this report. 
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TIME DETECTION AND MANUAL SORTING OF APPARENT STAR SIGHTINGS 


As shown in figure 4, the recorded coded star-mapper signals were compared with 
a threshold level and all signals exceeding this threshold were converted to a digital sig- 
nal with a star presence in a slit designated by a 1 and a star absence in a slit designated 
by a zero. Ideally, a nine-element code pattern of 110100010 would result from a passage 
of a single star across each group of slits. The last code element of 0 was added arbi- 
trarily to a single code group to provide a return to zero voltage condition at the termina- 
tion of each code group. The time duration T of a star signal generated by a star tran- 
sit across the narrowest transparent slit was set by 


w est 

where fi is 0.015° for each element of the nine-element code and the estimated spin rate 
w est was determined from examination of an oscillograph record of the raw star-mapper 
data in which time periods between repeated sightings of the more prominent star signals 
were assumed to be the time period for a complete vehicle spin cycle. 

The sample time used in the analog-to-digital converter was selected to be T/4 in 
microseconds. The presence of a star in a single element was indicated when two or more 
of four consecutive sample time intervals had signal levels which exceeded the variable 
threshold level. The uncertainty in the time of occurrence of a star presence or a star 
absence signal was T/2 microseconds. The output of the analog-to-digital converter 
was cross correlated in the decoder with a prestored replica of a code group as illustrated 
in figure 5. A star transit time t^ was acknowledged when correlation of the reference 
code sequence with any nine-element coded sequence produced agreement of at least two 
of four l’s and all five 0's or at least three of four l's and at least four of five 0's. The 
transit times of apparent stars were recorded T/4 microseconds after the time of occur- 
rence of the trailing edge of the ninth code element. The delay time of T/ 4 microsec- 
onds was necessary for readout of range time. The total uncertainty in star transit time 
included the ±5- microsecond uncertainty in range time and ±T/2-microsecond uncertainty 
in the time of occurrence of a star presence or star absence signal. 

At this point, a distinction must be made between the physical groups of slits on the 
actual reticle and the analytical slit edges for which star transit times are defined. A 
transit time, recorded by the correlator as the time at which a star transits some point 
fixed with respect to the reticle, corresponds to a line on the reticle located 9.25p after 
the leading edge of the vertical or slanted groups of slits. The analytical lines on the 
reticle described by 9.25 p after the leading edge of the vertical and slanted slit groups 
are defined to lie in planes designated as the vertical and slanted slit planes, respectively. 
A geometrical description of these slit planes is discussed in a later section. 
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The output of the correlator was a listing of apparent star transit times and their 
signal amplitudes. This listing was hand sorted in an attempt to relate transit times to 
discrete star sightings. The listing was examined for pairs of coded sequences of approx- 
imately equal amplitude such as a single star's transit of the field of view would be 
expected to produce. This amplitude sorting was repeated for all vehicle spin cycles 
as determined from the estimated vehicle spin rate. Additional sorting was done by 
searching for a recurrence of approximately equal time intervals between a pair of pulse 
sequences produced by the same star and the recurrence of approximately equal time 
intervals between pairs of sequences produced by different stars. The star-mapper data 
were then listed in the form of a pair of transit times for each apparent star sighting. 

STAR AND PARAMETER IDENTIFICATION 

The two primary operations of the star- mapper data processing procedure shown 
in figure 4 are the identification of the detected stars and the identification of the param- 
eters of the dynamical equations of motion describing the spacecraft attitude. The listing 
of star transit time pairs was operated upon in a star identification program which pro- 
duced a listing of identified stars and their sighting times. The identified star listing was 
used as input data to a parameter identification program, its resulting output being a set 
of the parameters which defines the dynamical equations of motion of the spacecraft. A 
priori data were used in both the star and parameter identification programs in the form 
of estimates of nominal vehicle trajectory, vehicle moments of inertia based on preflight 
measurements, and an approximation of vehicle spin rate based on successive sightings of 
identified stars. The following sections describe in detail the spacecraft attitude geome- 
try, the mathematical model assumed for the dynamical equations of motion, the star iden- 
tification program, and the parameter identification program. 

Space Attitude Geometry 

The general problem considered here is one of a geometrical formulation of the 
celestial attitude of a reference frame fixed in a spinning, torque-free, prolate, symmet- 
rical, rigid body. The problem is further defined as one for which the orientation of the 
body-fixed reference is to be defined at instants of time at which known stars are observed 
and as a function of time for periods during which no stars are observed. The basic mea- 
surements which are used to determine vehicle orientations are the times at which known 
stars appear in the transit slit planes of the star-mapper reticle. 

Basic reference frames.- The geometrical model for this problem involves the use 
of eight basic reference frames. Two additional reference frames were used to allow 
for star-mapper mounting misalinements and uncertainties in relationship of the vehicle 
principal axis to the vehicle body axes. The sequence of rotations and transformation 
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matrices for the defined Euler angles are listed in table I. The reference frames used in 
the solution are (see fig. 6): 

(1) Celestial reference frame with unit vectors iiJj^i- The unit vectors ii and 

ki are selected to be in the direction of the first point of Aries and the North celestial 
pole, respectively. The unit vector jj is selected to lie in the celestial equatorial plane 

to complete a right-hand system. 

(2) Inertial reference frame defined by two Euler angles $ and 0 relative to the 
celestial reference frame. In this frame the unit vector k 3 is defined as parallel to the 
angular momentum vector of the spacecraft. 

(3) A body-fixed reference frame defined by the Euler angles <p, 9, and ^ rela- 
tive to the inertial reference frame. In this frame the unit vectors ig> jg, and kg are 
defined to be parallel to the principal axes of the body. The unit vector kg is defined 
to be along the vehicle spin axis. Since the vehicle is assumed to be symmetrical, the 
unit vectors ig and jg may be selected to be any mutually orthogonal pair which is 
also orthogonal to kg. For convenience, the ig vector is considered to be directed 
along the nominal optical axis of the star mapper. 

(4) The reference frame of the star mapper relative to the spacecraft principal axes 
is defined by two misalinement Euler angles ei and e 2 - These two angles allow for 
misalinement of the star-mapper system relative to the spacecraft spin axis as measured 
prior to launch and also allow for changes in the location of true spin axis after launch 
and rocket fuel burnout. In this reference frame the ^ vector is selected to be along 
the optical axis of the star mapper. 

(5) The last reference frame of interest is defined by the location of the vertical and 
slanted slit planes (in which star transit times are detected) relative to the star-mapper 
optical reference frame. Two Euler angles y and /3 define this frame relative to the 
optical frame. There are actually two sets of y and 0; one set defines the vertical 
(y v ,0 v ) and the other, the slanted (y s ,/3 s ) slit planes. In either case the i 10 and k 10 
unit vectors are defined to be in the vertical or slanted slit plane. 

One additional angle is needed to complete the geometrical formulation of the prob- 
lem. This angle defines the coordinates of a star as sighted in a slit-plane reference 
frame. Figure 7 is a representation of a slit plane and the angles y and /3 relative 
to the optical axis of the star mapper. In this figure a star sighting is shown in the slit 
plane at an angle r) relative to the intersection of the ig - fs and i 10 ' k 10 P lanes - 
A unit vector to a star sighted in the slit-plane reference frame is 

r = cos rjiiQ + Oj iq + sin r} k 10 (!) 
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Figure 8 illustrates the geometry of the problem. For simplicity, the angles ej, € 2 , 

0y, Ty> rs> Vy> 311(1 Is are not shown - 

Preflight determination of j3 and y .- The optical axis of the star mapper was 
defined by a small transparent square called the alinement point, located in the center 
of the reticle. The relationship of this alinement point to the two groups of coded slits 
and to the radiometers and the azimuthal separation between the leading and trailing 
edges of both the vertical and the slanted slit groups were established by preflight mea- 
surements. These measurements were used to determine the four angles referred to as 
the slit parameter angles, /3y, /3g, yy, and yg. A discussion of the techniques used 
to obtain these preflight measurements is given in reference 3. A geometrical applica- 
tion of these measurements to the determination of the slit parameter angles is given in 
appendix A of this report. These angles are illustrated by the schematic of the reticle 
given in figure 9. The total angular relationship between the two groups of slits, defined 
by the angle /3, was determined. The angles /3y and /3g were defined as follows : 

/3y = 0 

As stated previously, the transit time pairs represent a known star's transit of 
points fixed with respect to the reticle. Schematically, these points lie on each of two 
analytical lines located on the reticle 9.25/1 after the leading edge of the vertical and 
slanted groups of slits. These two analytical lines were designated as the vertical slit 
plane and the slanted slit plane. The azimuthal angular relationships of the optical axis 
to these two analytical slits were defined by the angles yy and y g as shown in figure 9. 

Vehicle dynamics .- The relationship between the coordinates of an identified star in 
the celestial reference frame and the coordinate of the star as sighted in a slit-plane 
reference frame may be determined by use of the defined Euler angular rotations: 

R = (*)(©)(^)(0)(^)(e 1 )(e 2 )(y)O3)r (2) 

The vector R may also be written as 

R = cos a cos 6ij + sin a cos 6jj + sin 6kj (3) 

where a and 6 are the right ascension and declination of the identified star. The vec- 
tor r has previously been defined (eq. (1)) as 

A A A 

r = cos r/ijQ + °iio + sin ^10 
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where tj is defined as the elevation angle of the star in a slit plane. For the case of a 
spinning torque-free symmetrical rigid body, the angles <p and ip of equation ( 2 ) are 
defined as 

= <pt + 0 O W 

\p = \pt + \p Q (5) 


and the rates 0 and are defined as in reference 5. 


Ij cos 9 


zit ( 7 > 

: l 

where I 3 is the moment of inertia about the principal spin axis of the vehicle; Ij is 
the moment of inertia about either one of a pair of axes which are mutually orthogonal 
and orthogonal to the principal spin axis of the body; and w is the total spin velocity 
about the principal spin axis of the body. 

If all parameters of equation (2) were uniquely defined, the attitude of any set of 
axes would be known. However, only the y and /3 angles of equation (2) were assumed 
to be well known before flight. Of the remaining parameters in equation (2), ej, ^ e > 
rf> 0 , and ip Q were assumed to be unknown; estimates of # and 0 were available from 
launch data; and estimates of u>, <p, and 1 p were available from preliminary examina- 
tion of star-mapper data. With the use of identified stars, their sighting times, and least- 
square techniques, the parameters of equation (2) were identified. A summary of all a 
priori estimates based on preflight measurements and raw star-mapper data which were 
used in processing the data from Scanner flights 1 and 2 are listed in table II. 


Star Identification 

A star identification program was used to identify each star associated with each 
pair of transit times detected in the decoding routine. One of the primary requirements 
of the program was that the transit times be in chronological order and that the times 
appear in pairs, each pair being from the same star. The equatorial coordinates of the 
reference stars used in this program were obtained from the Boss catalog (ref. 6 ) and 
updated to the nearest 0.1 day of the time of flight by use of the astronomical techniques 
described in reference 7. The reference star listing was reduced by elimination from 
consideration all stars dimmer than a visual magnitude of +3.5 (ref. 3). The star listing 
was further reduced by selecting all stars falling within a selected band of the celestial 
sphere. The limits of this band, shown as an extended field of view in figure 10, were 
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defined to have the zenith angles of 76° and 104° away from the direction of spacecraft 
angular momentum vector as estimated from launch data. This band of ±14° about the 
spacecraft equator was necessary to accommodate the ± 3 ° field of view of the star 
mapper, an uncertainty of ±5° in the estimated angular momentum direction, and a pos- 
sible 6 ° vehicle cone angle. 

All stars of the reduced reference list were considered as candidates for the tran- 
siting stars. An additional listing of all these candidate stars was made in the form of a 
catalog of angular separation g c of each pair of stars in the reference list. Star identi- 
fication was principally a process of finding a cataloged angular separation which agreed 
to within a fixed tolerance Ag of an estimated angular separation g 0 of a pair of 
observed stars. 

A ngular separation of cataloged stars .- The angular separations of the cataloged 
stars are easily determined from their equatorial coordinates. In figure 11, the position 
of star G x on the celestial sphere is defined by coordinates (<*i, 6 i) and the position of 
star G2 is defined by the coordinates ( 02 , 62 )* The angular separation between Gj and 

G 2 , defined by arc GjG 2 and designated as g c , is given by law of cosines for spherical 
triangles as: 

cos g c = sin 6 j sin 62 + cos 61 cos 62 cos(aq - o^) ( 8 ) 

Angular separation of observed stars .- The angular separation between a pair of 
observed stars was estimated by use of the times at which the two stars transited the two 
slits of the reticle. This estimation was based on the assumption that vehicle cone angle 
and the misalinement angles were zero. These assumptions led directly to the definition 
of the azimuth-angle separation of a pair of stars as observed by the star mapper; that is, 

a = aJ est(tv(k+l) " tVk) (9) 

where w es t is the a priori estimated total spin rate and the subscripts k and V indi- 
cate the kth star transit of the vertical slit plane. The azimuthal separation of two 
observed stars as defined in equation (9) is analogous to the difference in the right ascen- 
sion angles of two cataloged stars as shown in figure 11. The elevation angle e k of a 
single star as observed by the star mapper may be written in a manner similar to that of 
equation (A4) in appendix A as 


tan e k = Sln & > est( t Sk ~ tyk) ~ (I Vyf + | rslfl , 

tan 0 S U 

where the subscripts Sk and Vk refer to the kth star transit of the slanted slit plane 
and the vertical slit plane, respectively. The estimated total angular separation g 0 of 
a pair of observed stars may be written in the identical form of equation ( 8 ) as 
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( 11 ) 


cos go = sin ek sin ek+i + cos ek cos ek+i cos a 

A pair of stars in the reference list was selected as the prime candidate for the identifi- 
cation of a pair of observed stars if 

I g c - go | = Ag 

where Ag is some fixed allowable difference in the two angular separations. This dif- 
ference was adjustable and was set to be that value which would result in the selection of 
a sufficient number of prime candidate star pairs so that at least three stars per vehicle 
spin period were identified. Positive identification of sighted stars was accomplished by 
identifying triplets of sighted stars by using the procedures discussed in appendix B. 

After star -identification processing was completed for one spin period, all informa- 
tion related to the identified stars was available. This information included the Boss 
catalog number, visual magnitude, right ascension, and declination of each identified star. 
The information also included a listing of transit times of each identified star for both the 
vertical and slanted slit planes. Each succeeding spin period was treated in a similar 
manner with no knowledge from previous spin periods assumed, since the flight data were 
expected to contain a minimum of three stars per spin period. After processing was com 
pleted for all spin periods, a synopsis made of all stars identified and the pertinent infor- 
mation for each star was printed. These data were next used in the parameter determi- 
nation program. 


Parameter Identification 

In the development of the spacecraft- attitude determination geometry, equation (2) 
was derived to show the relationship between the celestial coordinates of an identified star 
and the coordinates of the star in either the vertical or slanted slit plane reference 
frames. Substituting equations (1) and (3) into equation (2) results in 

/cos 6 cos a\ / cos ’A 

| cos 6 sin a j = ($)(0)((|))(0)(^)(ei)(e2)(r)(^)| 0 1 

\ sin 6 / \sln V / 

The parameters <p and i// have been previously defined (eqs. (4) and (5)) as 

4> = 4>t + <p 0 
\p = \p t + 

The nine undetermined parameters of equation (13) are 4>, 0, 4> 0 , 9, t//, \p 0 , ei, 

e 2 . An observed star will produce two equations such as equation (13), one each for the 
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slanted and vertical slit planes. If an adequate number of stars are observed, an iterative 
method may be used to solve for the undetermined parameters of equation (13). Since a 
large number of star observations were expected, a least-squares technique was selected 
for solution of the undetermined parameters. 

Equation (13) may be rewritten as r = JDEHR, where 

D =N" 1 ('0' 1 ' 

E = (^)- 1 ( 0)- 1 (^)- 1 I 

i i I (14) 

H = (0) _1 (<t>) -1 

J = (0 )' 1 (y)- 1 j 


The individual elements of the matrices J, D, E, and H are given in table m. The 

dependence of r on the undetermined parameters may be expressed in functional form 
as 

r(X)= JDEHR ( 15 ) 

where X is defined as the column vector of undetermined parameters 



The dependence of r on y and (3 is implied in equation (15), but these parameters 
were not included in the vector X since they were assumed to be well known before 
flight. The vector function r may be expanded about the true values of X by use of 
a Taylor series expansion as 


r(X) ~ l(X - AX) + F(X - AX)AX 


(17) 
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In equation (17) all terms higher than the first variational terms are assumed to be negli 
gible The matrix F is the matrix of first partial derivatives of the elements of r(Xj 
with respect to each element of X and evaluated at X - AX. Each element of F is 
expressed as 

9i 7 


F Zm = 


dx 


(18) 


m 


where 

and 


l - 1, 2, 3 
m = 1, 2, 3, . 


The vector X - AX represents an estimate of X, and AX represents the error in the 
estimate. This estimate of X is denoted by X hereafter. From the jth estimate for 

X, or X^,and k star transit times, equation ( 17 ) may be written as a set of k 

equations 

r(x) k * r(x^^) k + F(x ( ^) k AX (j) (19) 

The dimensions of the vectors and matrices in equation (19) are 

r (X) k - 3k x 1 

r(x^) - 3k x 1 


? (x (J) l 


(k — 1 , 2 , 3 , . . ■ , k) 


3k x m 


AX) - m x 1 


If the error in 
defined as Q k \ 


the approximation of the set of k equations represented by equation (19) is 


( 20 ) 


Qjf - r(x a) ) k - rfx^ + 

The quadratic form of the error is 

Q«> T qW - - £(X) k ] T [r(x 0> ) k - L(x) k ] ♦ [r(x«>) k - 

+ iX<j) T F(x 0) ) k j^(x ( ») k - r(X) k ] + ix(i) T F(x«>) k F(x0 ) ) k 4x(J> (21) 


r(j) 


16 


By using the usual method ot a least-squares solution, the quadratic form of the error Is 
“I"* “ S Par “ al " 6 WKh — *» “ - * e,u°a r , ‘to 

®fer Q k >TQ k ) - 0 - F (x«>)^(x<») k - r^] * f{x<J)) k F(x(J>) k ixG) (22) 

There are nine equations such as equation (22), one for each of the undetermined parame- 
ters. The error vector 4x0) may be determined from equation (22) as 


AX 


(j) _ 


F (xO))>(xO))T^(xO))T E (x<i)) . 4 x) 


(23) 


A new estimate of X may now be determined 


as 


x (i+1 > = x<i> + axO) (24) 

»here X<J) > is the jth estimate of X and AX<i> is the solution to equation (23) This 

sta~ S "it ““ ° f ‘ he VarlaWe “ a * - -ch slit for eac 

star transit. Since the elevation angle was variable and not measurable, an alternate for 

mulation of the problem was desirable for simplicity of problem solution. 

star JanluTl 0 " 1“ ' 13> Sh ° WS lha ‘ ‘ here 3re tW ° dependent equations for each 
inden. T ^served star results In a pair of transit times and a pair of 

pendent equations for the slanted and vertical slit plane transits. If the parameters 

rr° n r were known ’ ihe equa “° as £ ° r the sianted sat ^ ^it ^ P re. 
sr rri, ° f the e<!uations resu,ung fr ° m tte ver “ cai siu ^ 

time This condition means that knowledge of „ Implies knowledge of „ and this 

condition further implies that only two independent equations exist for each pair of star 
transit times. The problem solution can therefore be simplified by making use of the 

constraint equation represented by the second element of equation (13). Applying this 
constraint to equation (15) results in 


0 = J2DEHR 


(25) 


where the subscript 2 on J 
previously described method to 
vector 


refers to the second row of the matrix J. Applying the 
equation (25) results in a new definition of the correction 
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where the subscript 2 refers to the second row of matrix F and the second element of 
the vector r. Dimensions of the matrices in equation (26) are 

AX^) — m x 1 


F 2 (x (j) ) k - k x m 

^X - k * 1 


As before (eq. (24)), a new estimate for X may be obtained by 

xti+D = xf j) + Axf^ 


where X (j) is the old estimate and AX^ is the solution to equation (26). The updating 
of the estimates of X may be continued until each element of AX<1> S M, where M s 
a predetermined convergence criteria based on the required accuracy of convergence. 

If either or both of the parameters 9 and 0 are zero, the matrix 


F 2 (xa)) k ] T [F 2 (xO)) k ] 

of equation (26) is singular and no solution for AX is possible. If these two parameters 
are nonzero, but nearly zero, the matrix Is nearly singular and special care must be 
taken to invert this matrix. In the case of the Scanner vehicle launchings, the angle © 
was never less than 50° and therefore presented no problem. However, the vehic e cone 
angle 6 was expected to be less than 6° and special treatment was necessary for inver- 
sion of the matrix. A discussion of the inversion procedure Is presented in appendix C. 


Initial Estimates of X 

Before a solution of the values of the nine undetermined parameters of X is possi- 
ble, initial estimates are needed. These initial estimates were obtained by considering a 

reduced problem of the three undetermined parameters i // 0 + 4> 0 , and ©• 

assuming 0, € 1 , and e 2 to be zero, <f> and * were calculated by use of w e st and 
the preflight measured values of Ii and I 3 . The values 01 u an » 
this solution were estimated from launch data. The initial estimate of <p Q + ^ 
three-unknown problem was determined by using equation (13) with the assumption tha 
yv was zero, and setting the vertical slit transit time of the first identified s ar o e 
Uro. m addition, the elevation angle was assumed to be zero. Equation (13) then 

becomes 
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'cos 6 cos 


1 


cos 5 sin a 1 = ($)( 0 )(<£ o )(i// o ) 0 
sin 6 / W 


( 27 ) 


By making use of the definitions of the transformation matrices listed in table I, equa- 
tion (27) becomes 


/bos 6 i cos aq\ /cos $ -cos 0 sin <f> sin 0 sin Jcos(<p 0 + \p 0 f\ 

cos sin cki = sin $ cos © cos $ -sin 0 cos $ sin( 0 o + \ p 0 ) 


(28) 


sin 6 i 


sin 0 


cos 0 


where the subscript on a and 6 denotes the right ascension and declination values of 
the first observed star. The matrix equation (28) may be solved for cp 0 + i // 0 as 


tan (4>q + i/'q) 


-cos 6 i cos 0 sin(<l> - aq) + sin 0 sin 6 j 
cos cos(4> - aq) 


(29) 


In this three-unknown problem, the AY correction vector as defined in appendix C 
becomes 



where the elements Ay 2 , Ay 4 , and Ay 5 were computed by the least-squares program 
discussed previously. The solution for X resulting from this procedure produced an 
initial estimate for \}/ Q + cj) 0 and better estimates for $ and 0 which were then 
applied to the problem of the nine undetermined parameters. This approach produced 
satisfactory estimates of the sum of \]s 0 and <f> 0 , but no knowledge of estimates of their 

individual values. The procedure used here was arbitrarily to set 4> 0 = ^ 0 = 0 

as initial estimates for the nine-unknown problem. This procedure is a reasonable one 
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since the solution of the nine-unknown problem is heavily dependent on <p 0 + \p 0 and 
lightly dependent on their individual values for a small cone angle as discussed in appen- 
dix C. A summary of all a priori estimates is listed in table n. These a priori estimates 
were used to initiate the three-unknown problem. 

DEVELOPMENT OF DATA-REDUCTION ACCURACY CRITERIA 
FROM PRE FLIGHT ANALYSIS 

Before application of the attitude determination technique to flight-data reduction 
the accuracy of the parameter identification method had to be ascertained in order to 
establish expected accuracies for the Project Scanner flights. In addition, criteria had 
to be developed in order to assess the adequacy of any estimate of X as determined 
from flight data. To establish the expected accuracy of spacecraft attitude determination, 
a simulation study was performed. Criteria for assessing the adequacy of estimates of 
X were developed by using the results of this study. In this study simulated flight data 
were generated by using the expected launch and flight conditions of Project Scanner 
flight 1. The mathematical model of simulated spacecraft equations of motion was iden- 
tical to that which was assumed for the flight vehicle. Simulated flight data were in the 
form of a time list of transited stars and the right ascension and declination values of 
these stars updated to the expected launch time. These simulated data were used to esti- 
mate the values of the previously described parameters of the vector X in the same 
manner as for flight data. The parameters of X were used to generate pointing- direction 
time histories of selected axes. The errors in these pointing directions, specifically the 
pointing-direction errors of the optical axis, were the basis of assessing the accuracy of 
the attitude determination solution. 

Pointing-Direction Analysis 

Figure 12 shows a flow diagram illustrating the generation of time histories used in 
the simulation studies. A set of nine parameters defining the equations of motion, launch 
data, the reticle configuration of Scanner flight 1, and the right ascension and declination 
of cataloged stars were used to generate simulated star transit-time data. The simulated 
transit times plus expected errors and the right ascension and declination values of the 
transited stars were processed by the parameter identification program in the same man- 
ner as flight data. In addition, the a priori estimates used in this processing were identi- 
cal to those used for flight data processing. The values of the nine parameters deter- 
mined by the parameter identification program along with the known values of 0y, /3g, 
yy, and yg were used to generate time histories of the pointing directions of the space- 
craft principal spin and the star-mapper optical axes. The variation of error in these 
pointing directions as a function of several parameters was investigated and is presented 
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in the following sections. These parameters include the length of time covered by a least- 
squares data fit in terms of fractional parts of a precession cycle, number of stars per 
spin period used in the data fit, and angular separation of stars within a spin period. 

Generation of error time histories.- Since the information provided by the star 
mapper was transit times of observed stars, the basic errors in the system were inac- 
curacies in the detection of transit times for the vertical and slanted slit planes The 
two sources of error which contributed to these transit time errors were the inaccuracies 
in time detection of a star transit across the slit planes and resolution of system clock 
time. As previously discussed, these two error contributions were ±T/2 microseconds 
and ±5 microseconds. The total time error was conservatively assumed to be the sum 
of the two errors or ±1 + 5 microseconds. This total time error was assumed to have 
a uniform probability density function and therefore a standard deviation of approximately 
19 microseconds for the nominal value of T being approximately 56 microseconds. In 
the simulation studies, errors having a standard deviation of approximately this value were 
added to the transit-time data, by selecting, with equal probability, error values of 26 0 
or -26 microseconds. ' ’ 


Different sequences of random time errors were added to each selected sample size 
of the simulated data. Thus, for each sample size of star data, a run through the least- 
squares program was necessary for each noise sequence to determine the resulting esti- 
mates of the original nine parameters specifying vehicle dynamics. Each set of nine 
parameters was then fed into a time-history generator, that is, the vehicle dynamic equa- 
tions written as a function of time. Beginning at time zero, and for equal intervals of 
time, a time history of attitude data in celestial coordinates (right ascension and declina- 
tion) for the vehicle symmetry axis k 6 and the star-mapper optical axis i 8 was calcu- 
lated for each set of nine unknowns. Each of these time histories was computed for the 
length of time covered by the selected sample size of simulated data used in the least- 
squares program. This computation provided time-history records of attitude both at the 
times of star sightings and during the periods of time between star sightings. The atti- 
tude data representing times between star sightings provide a measure of the interpolation 
capability of the attitude-determination procedure. Each time history associated with a 
transit-time noise sequence was compared with a reference time history generated by 
use of the original nine parameters used in the simulation program. This comparison 
resulted in an ensemble consisting of as many error time histories in right ascension and 
declination for the k 6 and i 8 axes as the number of noise sequences used in the genera- 
tion of pointing- direction time histories. ° 

Pointing-direction errors.- Standard-deviation calculations of the errors in right 
ascension a a (t v ) . and declination a 6 ( t v ) for each of these axes were made across the 
ensemble at specific instants of time t v . For example, a right-ascension standard- 
deviation calculation at a specific instant of time is designated as cr^t^g and a^t^g 
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for the vehicle principal axis and the star-mapper optical axis, respectively. The 
expected standard- deviation error in pointing direction cr a (t„) of the k 6 or 8 «« s 

was assumed to be the root sum square of the standard deviation of the errors in rig 
ascension and declination of each of these axes. The maximum value of this root sum 
square value calculation is defined in this report as o a (t„) max and was determine y 

computing the value of a a (t„) at each t„ and then selecting the maximum value of 

a (t ) Spacing of the time points t„ was examined for 0.25 second and 0.01 second, 
negligible differences being noted in the standard -deviation calculations. The root-mean- 
square value of the pointing-direction error was calculated for the time duration of sev- 
eral error time histories. It is interesting to note that the results of these calcula on 
agreed with standard- deviation calculations made on an ensemble basis to within approxi- 
mately 15 percent. The use of ten sequences of data for these calculations and applica- 
tion of the chi-squared test resulted in a 99-percent confidence that the true value o 

CT a (ty) m ax would lie ln the range 

0.618CT a (tj,) max = true CT aM ma x = 2 - 278a a(^)max 

g.„^„rv of pointlng-direction-error analysis .- The values of the nine parameters 
used in the generation of the simulation data were 

0 =0.310758° 


xP o = 51.081051° 
ip = 287.844975°/sec 
d> = 316.572701° 

0 = 19.137575°/sec 
$ = 76.462935° 

© = 54.126671° 


€l = 0.064170° 
e 2 = 0.031017° 

The vehicle inertia values were assumed to be 

I x = I 2 = 129.92 slug-ft 2 (176.15 kg-m 2 ) 

I 3 = 8.27 slug-ft 2 (11.2 kg-m 2 ) 
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A partial listing of the simulated star transit data is 


SIMULATION DATA LIST 


Star 

Transit time, sec 

Visual magnitude 

Boss catalog number 

Vertical slit 

Slanted slit 

154 

3.908927 

3.912402 

3.000 

19607 

172 

3.966026 

3.971817 

2.310 

20947 

177 

3.979605 

3.996827 

2.000 

21491 

219 

4.183629 

4.193429 

2.940 

25180 

221 

4.196009 

4.211718 

3.300 

25661 

225 

4.211204 

4.229113 

2.710 

26161 

239 

4.293068 

4.304024 

3.210 

28682 

23 

4.503505 

4.505614 

3.420 

3584 

43 

4.652463 

4.653949 

2.920 

6274 

48 

4.690902 

4.702167 

1.700 

6668 


The star number is the order of occurrence for the star in the Boss catalog. 

Information from this listing which is used by the least-squares program includes 
the star index number and the transit times for the star across the two slits. The star 
number was used for a storage index number of the right ascension and declination values 
of the particular star transited. 

The accuracy of the attitude -determination procedure was determined by calculating 
°a^^)max * rom data representing multiples of a precession period, number of star 
sightings per spin period, and angular displacement of stars within a spin period. The 
results of these calculations showed that star placement within a spin period was a major 
factor, since with poorly placed stars, for example, all stars in a spin period are placed 
within half a period, precession information is not well defined. As indicated in figure 13, 
the accuracy of attitude determination increased as the data-fit period used in the least- 
squares program increased and as the number of tar sightings per spin period increased 
for the one-precession- cycle case. The one-half precession- cycle case also indicates 
evidence of this trend. The one-and-one-half precession-cycle case shows no particular 
trend and the results for this case were interpreted to mean that the errors in the esti- 
mates of X were in the noise level for any considered number of stars per spin period. 
Based on the plots of figure 13, the predicted Scanner star mapper accuracy was esti- 
mated to have a ff a (tj,) max value of approximately 0.004° since most of the least - 
squares-data fitting covered flight-data periods of three-fourths of a precession cycle 
and contained at least 4 to 6 star sightings per spin period. The values of °a(Mmax 
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obtained in the simulation study gave assurance that the attitude-determination technique 
would meet the experimental requirement of a <r a of 0.016° with a confidence level of 
99 percent. 


Relation of Transit-Time Errors to Least-Squares Functional Residuals 


To evaluate the accuracy of attitude determination resulting from use of the esti- 
mates of X as computed in the least-squares program, time errors e ^ were added to 
the simulated star transit times tfc. The values of etk had an equal probability of 
having values of +26, 0, or -26 microseconds and a standard deviation a t of approxi- 
mately 21 microseconds. Since it is not possible to measure for flight data, it was 
necessary that some other parameter be used to estimate the magnitude of transit time 
errors. The functional residuals of the least-squares solutions resulting from simulated 
data were examined for this purpose. 


The adequacy of any estimate of X may be determined by the use of certain 
assumptions. If is assumed to be a good enough estimate to allow the use of 

xO+l) and true kth transit time tfc. to satisfy approximately the constraint of equa- 
tion (25), then 

J 2 DEHR ~ 0 (31) 

For convenience, let the functional dependence of equation (31) on X and t for any star 
and the (j + l)st estimate be denoted by 


J 2 DEHR = f(x^ +1 \ 4) 


(32) 


Equation 

4 as 


(32) may be approximated by a Taylor's series expansion about the true time 



(3+D t t 

’ l k 



(33) 


By defining 

4 ~ At k = tfc ( 34 ) 

as any observed transit time, and by setting f(x(j +1 ),4) equal to zero by use of equa- 
tion (31), the time error in the observed kth star transit time may be approximated by 


Atk 



(35) 
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The convergence factor M was adjusted until a least-squares iteration resulted in 
time residual values At^ which agreed with to within an accuracy of 5 microsec- 
onds. In this study the standard deviation of Atjj was calculated for each set of k 
stars used in a least-squares solution. This standard deviation of the time residuals 
qt r was found to be approximately 20 microseconds for all cases in the study. For 
cr^ r ~ 20 microseconds, the accuracy of the computed vehicle-attitude orientation was well 
within the specified requirements. This parameter a^ r was the factor which was used 
to assess the adequacy of a least- square solution for flight data. In the case of flight 
data, the convergence factor M was adjusted until aj r was near its minimum value. 
When ofc r was at an acceptable minimum value, a least-squares solution was accepted 
as an adequate one. 

Individual values of Atk were also used to assess the adequacy of each least- 
squares solution using flight data. If one or more values of Atjj in a set of k data 
points were inordinately larger than the remaining members of the set, the questionable 
data points were reexamined. In every case the large values of At^ could be attributed 
to errors in range time counting, incorrectly cataloged right ascension or declination 
values of the star associated with the particular data point, or false star identification 
resulting from two or more detectable stars appearing in the star- mapper field of view 
simultaneously. The errors in range time were relatively easy to locate since a taped 
time-history record of range time counting was available. The errors in cataloged star 
data were determined by cross checking several astronomical tables of star listings. The 
identity of the incorrectly cataloged stars and the procedure for handling multiple stars 
in the field of view are discussed in the section "Flight Results." 

Extrapolation of Attitude Determination 

The star-mapper data-processing procedure was further studied to determine the 
feasibility of predicting or extrapolating in time the pointing directions of the ig and 
kg axes based on the results of the least-squares program. In this study the extrapola- 
tion period covered a period of time approximately 1 precession cycle after the period of 
time used in the least-square solution of the nine parameters defining the vehicle equa- 
tions of motion. The ability to predict or extrapolate attitude information was limited by 
transit time errors and the convergence criterion used in the least- squares program. 

The two principal members of the set of estimated nine parameters which affected the 
ability to predict attitude were 0 and 0. The region of extrapolation was also governed 
by the allowable magnitude of cr a (tj,). Inconclusive results were obtained in this study and 
extrapolation for more than a fractional part of a precession cycle was not possible. It 
was therefore decided not to attempt attitude extrapolation by use of Scanner flight data. 
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FLIGHT RESULTS 


The star mapper was flown on two suborbital flights of the Scanner vehicles to pro- 
vide data necessary for determining celestial attitude of these vehicles. On both flights, 
a cold-gas-reaction jet-control system was used to erect the vehicle to within ±2° of the 
local vertical at four definite times during the flight in order to maximize the data gath- 
ering of the primary-horizon-definition experiment. During the intervals between these 
erection times, the controls were turned off and the spin- stabilized vehicle was assumed 
to be free of external and control torques. Star mapper data for each of the two Scanner 
flights consisted of four such data periods, each of which was approximately 2^ minutes in 
duration. Had the vehicle actually been a symmetrical, rigid, torque-free, spinning body, 
the vehicle's motion for an entire data period could have been represented by a single set 
of nine vehicle motion parameters. Under this assumption the single set of nine vehicle 
motion parameters could have been determined by performing a least-squares fitting of 
all the star transit times in an entire data period. Preliminary evaluation of the results 
of the attitude determination procedure for both flights indicated that some undefined 
source produced small disturbing torques during the data periods. 

Desired accuracies in the determination of vehicle attitude were obtained with the 
assumed model by reducing the time periods covered by a least- squares solution to 
approximately three-fourths of a precession cycle in order to minimize the effects of 
these torques and possible vehicle asymmetries. For example, the second data period 
of Scanner flight 1 was divided into 12 time intervals and 12 sets of nine vehicle motion 
parameters were used to define the vehicle's motion during that data period as shown in 
table IV. The discrepancy between the model assumed for the motion of the vehicle in 
the parameter-identification program described earlier and the actual vehicle motions 
may be seen by observing the total angular spin rate o>, which is defined by 

co = + (p cos 9 

For a symmetrical, rigid, spinning, torque-free body, this angular spin rate should remain 
constant. Spin rate determined by succeeding solutions of approximately three-fourths of 
a precession cycle over an entire data period was slowly but constantly increasing. This 
constantly increasing spin rate was characteristic of all four data periods of both Scanner 
flights. The spin rate obtained for the second data period of flight 1 is shown in figure 14. 
The increase of co over the time period of the data period of figure 14 was 0.04 percent. 
This percentage of increase for co was approximately the same for all data periods of 
both flights. No attempt was made to determine the source of this variation in the 
assumed model. A possible source of the small disturbing torques could have been 
incomplete closure of the control valves in the reaction-jet control system during the 
data period. 
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The number of stars per spin period sighted by the star mapper during each data 
period varied from 3 to 14. The listing of identified stars sighted by the star mapper 
technique during the Scanner flight 1 is shown in table V. Ambiguities in associating raw 
data with vertical or slanted slit transits and in star identification were encountered when 
two or more stars with small angular separations were present in the field of view of the 
star mapper. Such problems usually resulted in unusable data from one or both stars, in 
which case the questionable stars were eliminated from the star identification list. As a 
result, the number of stars actually used for attitude determination of the Scanner vehi- 
cles, shown in table VI, was considerably smaller than the number identified. 

Preliminary data results revealed that several stars consistently exhibited large 
time residuals. Subsequent checking traced the cause of these high residuals to errors 
in the updated star positions based on values listed in the Boss catalog. The positions 
of these stars (6 Orionis, rj Ursae Majoris, and a Piscis Austrini) were cross checked 
with positions taken from the apparent place tables (ref. 8) and updated to the nearest 
0.1 day of the flight. Discrepancies of the order of 0.01° in the star positions taken from 
the Boss catalog and updated to the nearest 0.1 day of flight were found. When the stars' 
positions determined from the apparent place tables were used for these stars, the time 
residuals for these stars dropped to acceptable limits. 

An analysis of all solutions of the nine vehicle-motion parameters determined from 
flight data showed a time residual c^ r of approximately 30 microseconds. In the simu- 
lation study described in an earlier section, star-transit-time errors having a standard 
deviation of approximately 21 microseconds were added to the simulated star transit data. 
These time errors resulted in a standard deviation of the time residuals cr^ r of approxi- 
mately 20 microseconds when data were fitted over an interval of 0.5 to 1.50 precession 
cycles. The corresponding maximum error in the pointing direction cr a was found to 
be approximately 0.004°. It can be shown analytically that there is a linear relationship 
between o a and cr t . Since it has been shown that oj. and oj. r have very close agree- 
ment, it can be said that there is approximately a linear relationship between a a and c^ r . 
Therefore, the expected accuracy in pointing direction determined from flight data is 
approximately 0.006°. Although the time residuals were somewhat larger than the tran- 
sit time errors used in the simulation study, the accuracy of the data still remained within 
acceptable limits for the flight experiment. 

To illustrate the use of the nine vehicle- motion parameters, a set of these parame- 
ters was selected from the second data period of Scanner flight 1 and was used to generate 
time histories of the star- mapper optical axis ig and the vehicle principal spin axis kg. 
The results of these computations are shown in figures 15 and 16 in which the pointing 
directions of these two axes are plotted in terms of declination and right ascension for 
selected time periods. Figure 15 represents the pointing direction of the optical axis 
during one spin period and figure 16 represents the pointing direction of the vehicle 
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principal spin axis during a precession cycle. Figure 17 represents the pointing direction 
of the kg axis during a precession cycle. The vehicle motion parameters used to gen- 
erate this plot were the same ones as used for figure 16. Figure 17 illustrates the effects 
of the misalinement angles, and 63 * The complete listing of all the sets of nine 
parameters for each data period and for both flights was used in a similar manner to 
generate time histories for the i 8 axis. These results along with the preflight mea- 
sured angular relationship between the star-mapper ig axis and each radiometer opti- 
cal axis were used to determine the time history of each radiometer optical axis in celes- 
tial coordinates. 

POSSIBLE EXTENSIONS OF THE METHOD OF 
ATTITUDE DETERMINATION 

If the reported attitude-determination technique is used in future space flights, sev- 
eral areas of the data-processing procedures should be modified to facilitate data han- 
dling. The extent of the modifications depends on the particular application. However, in 
any future application, it would be desirable to eliminate, as much as possible, manual 
data handling, particularly in the areas of pairing of apparent star observation times and 
star identification. Increased automation in these areas may require that the entire data- 
processing procedure be a closed-loop operation. For application to an orbiting space- 
craft, a more exact model of the dynamical equations of motion would be required. This 
model should include the effects of vehicle asymmetries and at least the most significant 
external disturbance torques for the spacecraft under consideration. Methods other than 
the reported least- squares method used for the solution of the undetermined parameters 
in the equations of motion should be investigated. In particular, sequential updating meth- 
ods such as the Kalman filter should be studied for application to long-life orbiting 
spacecraft. 

A least-squares solution may be improved by using the time residuals to correct for 
errors in star transit times. Use of these time residuals to correct star-transit-time 
data is expected to result in improved least-squares solutions and therefore improved 
attitude- determination accuracy provided that the model of the vehicle equations of motion 
is well known. If greater attitude determination accuracies were required for the Scanner 
flights, this method could be investigated for possible improvement of the achieved accu- 
racies. The value of this method would first need to be evaluated for the cases in which 
the true vehicle dynamics and the assumed model dynamics differed. 

If the transit time errors are truly random, the time residuals will exhibit no pre- 
ferred periodicities. In the case of the Scanner flights, periodicities which agreed with 
the periodicities of the radiometer scan mirrors were observed in the time residuals. 
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This result was expected since preflight analysis of the effects of motion of the radiometer 
scan mirrors indicated that the vehicle moments of inertia would vary with a period equal 
to the period of the scanning mirrors. The maximum predicted change in the moments of 
inertia based on the mirror scan motion was 0.017 percent. Simulation studies indicated 
that asymmetries of this amount would produce negligible error in attitude determination 
when a symmetrical model was assumed. The important point is that the time residuals 
of the Scanner flight contained periodicities which were directly attributable to ignored 
but known sources which contributed to the vehicle dynamics. Thus the time residuals 
might be used to reveal periodicities attributable to unknown sources. This information 
could then be used to define a more exact dynamical model. 

CONCLUDING REMARKS 

An attitude-determination technique using star-transit information has been suc- 
cessfully demonstrated with an estimated accuracy greater than 0.01° during suborbital 
flights of the spin-stabilized Scanner spacecraft. Identification of the sensed stars was 
made on the basis of their angular separations in the celestial sphere as determined from 
their transit times. Information consisting of the cataloged positions of these identified 
stars and their measured transit times was fitted in a least-squares method to a set of 
equations representing an assumed dynamic model of the spacecraft attitude. The 
dynamical model was developed on the assumptions that the vehicle was spin-stabilized, 
symmetrical, and free of external torques. The least-squares fitting of star position 
and transit time data resulted in definition of sets of parameters of the assumed dynam- 
ical model. The equations of the model and the sets of parameters were then used to 
generate time histories of selected reference frames associated with the model. In par- 
ticular, time histories were determined for pointing directions of the vehicle principal 
spin axis and the star- mapper optical axis. Accuracy of the solution of the attitude deter- 
mination problem was based on the pointing-direction accuracy of the star-mapper optical 
axis which was the only axis that could be related to the radiometer lines of sight. 

Estimated accuracies of the pointing- direction determination of the star- mapper 
optical axis from flight results agreed favorably with predicted results. The achieved 
accuracies were considered to surpass the results obtained from other attempts to mea- 
sure or determine celestial attitude of a spin- stabilized spacecraft. 

Differences between the predicted results and actual results were evidenced by 
greater transit-time residuals in the flight data. These larger time residuals may be 
attributed to either an increase of noise or to an inexact model of the assumed vehicle 
dynamics. The latter reason is considered to be the most probable since it was observed 
that total spin rate steadily increased as a function of flight time. A source of this 
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problem might have been incomplete closure of the control valve in the reaction- jet con- 
trol system. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., May 6, 1968, 

715-02-00-01-23. 
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PREFLIGHT MEASUREMENTS OF THE RETICLE AND DESCRIPTION 
OF THE DETERMINATION OF THE SLIT PARAMETER ANGLES 

The angular relationship between the slit plane, defined by the ijo " ^10 plane, 
and the ijQ - kg plane of figure 7 is defined as the angle /3. The actual experimental 
measurements of the reticle configuration have been described in reference 3, and no 
attempt is made to redescribe these measurements. However, since an understanding of 
the geometry of this measurement is essential to any discussion of the method of deter- 
mining the angle ( 3 , a summary of the method described in reference 3 is discussed here. 

A collimator beam was used as a light source and the star mapper was mounted on 
a table which could be moved in azimuth and elevation. The table could be locked in ele- 
vation and moved in the azimuth direction. As nearly as possible, the vertical slit group 
was fixed so that it was perpendicular to the azimuthal plane of the table. The nonverti- 
calness of the vertical slits to the table azimuthal plane is accounted for in the computa- 
tion of /3. The star mapper and table were positioned such that the collimator beam was 
at the center of the alinement point. The alinement point defines 0° in the elevation and 
0° in the azimuthal plane of the table. 

With the table locked at an elevation of 0°, the table was rotated in the azimuthal 
plane until the collimator beam was coincident with the leading edge of the vertical slit 
group at a point defined as V Q . The azimuthal angle of this rotation, defined as p Vo , 
was recorded. 

The table was then returned to the alinement point (0°,0°) and rotated in the oppo- 
site direction until the collimator beam was coincident with the leading edge of the slanted 
slit as S 0 and the azimuthal separation between the alinement point and the leading edge 
of the slanted slit, defined as p So , was recorded also. In figure 18 ON defines the axis 
about which azimuth was measured and P is the alinement point which defines 0° eleva- 
tion and 0° azimuth. 

When the table was returned to the (0°,0°) position and elevated to an angle 4, the 
elevation axis was then positioned at with NP^ defined by an arc of length (90° - 4) 
as shown in figure 18. The azimuthal angle p g ^ was measured by the rotation about 
ON from P^ to St. 

The spherical angle S^NS 0 is denoted by A where 

x = Ps£ " PSo 
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The angle /3j_ defines the angle of nonverticalness of the slanted slit group with the 
azimuthal plane of the table. 

By the law of sines for spherical triangles, 


cos £ sin A 

sin /3j = 7 =^ — 

sin S 0 S| 

Since the arc S^N equals 90°, the law of cosines reduces to 


COS 01 = 


sin £ 


sin S^Sq 

Dividing equation (A2) by equation (A3) reduces to the form 

sin A _ s * n (^S£ PSo) 


tan 0 i = 


tan £ 


tan £ 


(A2) 


(A3) 


(A4) 


This same process was repeated for the vertical slit group to obtain a corresponding 
angle 02 between its leading edge and the plane perpendicular to the azimuthal plane 
of the table at point V Q on the leading edge of the vertical slit group. The resulting 
equation, similar to that for 0 ^, defines 02 


tan 02 = 


s^ (Py£ - Pyp) 
tan £ 


(A 5) 


The angle between the slanted slit plane and the i^Q - kg plane of figure 7 is defined by 
the angle 0 g as 

0 S = 01 + 02 (A6) 

and the angle between the vertical slit plane and the iiQ - kg plane is defined as 

0y = 0 (A7) 

In each case the slit plane is defined as the iiQ - kig plane of figure 7. Preflight mea- 
surements were also made to determine the azimuthal separation between the leading and 
trailing edges of the vertical and slanted slit groups. 

The azimuthal angular relationship of the vertical slit plane to the optical axis is 
defined as 

Yy = Py o " (A 8 ) 

The azimuthal angular relationship between the slanted slit plane and the optical axis is 
defined as 

Y s = P So ~ 9 - 25 M (A9) 
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The preflight measurements indicated a different value of p for the vertical and slanted 
slit groups. The measured values of p for the vertical and slanted slit groups were 
used in equations (A8) and (A9), respectively. 
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IDENTIFICATION OF STAR TRIPLETS 

Recall from figure 4 that the star-mapper data at the input to the star-identification 
program were in the form of star-transit-time pairs where each transit-time pair 
resulted from observation of a single star. This time listing was grouped in triplets of 
transit-time pairs for each vehicle spin period which was obtained from u> es t. This 
grouping is illustrated in the following example sketch: 

Triplet 1 


Star 1 Star 2 Star 3 Star 4 Star 5 



tyi tsi *V2 ts2 tv3 ts3 *V4 ts4 tv5 ts5 


Spin period 

The first triplet is composed of transit-time pairs 1, 2, and 3 and the second triplet is 
composed of star-transit-time pairs 2, 3, and 4. This grouping into triplets was contin- 
ued until all transit-time pairs in a spin period were exhausted. 

The angular separation of each pair of stars in a triplet was determined by equa- 
tion (11). These angular separations were used to obtain prime candidates for the 
observed star pairs from the reference list of stars. As many sets of prime candidates 
were selected for stars 1 and 2, then 1 and 3, and finally 2 and 3 as could satisfy equa- 
tion (12). If this is done properly, at least one of the prime candidates for star 1 of (1,2) 
and star 1 of (1,3) should be the same. This procedure was repeated to identify star 3 of 
(1,3) and (2,3), and star 2 of (1,2) and (2,3). 

This procedure is illustrated by the following example: Given a computed angular 
separation of two observed stars, the listing of angular separations of all pairs of stars 
in the reference list was searched for angular separations meeting the condition of equa- 
tion (12). Each pair of observed stars was associated with a sublisting of prime candi- 
dates taken from the reference star list. After this procedure was completed for 
stars (1,2), (2,3), and (3,1), the following sublists of prime candidates were available: 
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he star numbers denoting candidate stars merely refer to the index numbers chosen for 
he reference list of stars and should not be confused with Boss catalog numbers. Simul- 
taneous comparison of the three lists of prime candidates should yield one triplet of ref 
erence stars for the triplet of observed stars 1. 2, and 3. In the example, it ^ apparent 

7; h r y 7 ° f reference stars whlch has ““">011 stars in each of the three lists 
of candidates star pairs is 7, 9, and 29. This comparison results in a tentative Identifl- 

lon o o served stars 1, 2, and 3 as reference stars 7, 9, and 29, respectively. 

trinlet S o, m reJ meS * *?“ °' ° bSerVed StarS WaS identifiable hy more than one 

each e r n r ' The amMgUlty was overcome by using the identified stars of 

each possible triplet and the associated six transit times in a least-squares program to 

solve for the reduced set of parameters * + 0, and *. All triplets were rejected 

“r° f f ^ * d l sa greed by more than ±10° from the launch-data estimates 
o these two angles. The triplet which produced the smallest functional residual in the 

lTca^dTd a a"s. Pr0gram W3S Sel6Cled “ ide ” tified ‘ riP,e * ° U ‘ °' *“ lhe tr ‘P- 


Once a triplet was tentatively identified, the index numbers of the identified refer- 
ence stars were stored in a sample table as follows: 
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Identification of the second triplet stars ideally reproduces the index numbers 1 and 
16 since the next triplet is obtained by dropping the first star and adding a new one at the 
end. In the example, the added star is assumed to be index number 38 which yie s 



If the program was unable to identify the stars in a triplet, a negative number ( 1, 
-2 or -3) was stored in the table in place of a reference star index number, as indicated 
in the following list representing the identification of all stars in a complete spin period: 



The final step in the identification process was to examine the reference star index num- 
bers in all but the last two columns of the table of identified stars for identical index 
numbers. If two or more index numbers in each column were identical, then positive 
identification was assumed for the reference star associated with that index number. If 
no two index numbers in a column were alike, identification was not made. Identification 
of the first and second star in a spin period required that the first two transit-time pairs 
be updated by a spin period, the star pair representing these modified transit times be 
repeated at the end of the list, the second and third entries of the next to last column e 
copied into the first column, and the third entry of the last column be copied in the sec- 

ond column. 

The authors wish to acknowledge the contribution of Bruce Vannelli of Control Data 
Corporation who, in the course of writing a digital computer program for star identifies- 
tion, developed the logic used in the identification of stars. 


♦ 


36 




APPENDIX C 


LEAST- SQUARES FORMULATION FOR SMALL CONE ANGLES 


Equation (25) may be expanded in a Taylor series expansion as 

0 ~ J2 DEHR + F 2(x)aX (Cl) 


with terms higher than the first variational terms assumed to be negligible. The matrix 
F 2 (x) is defined as 


Fofx ) = (j- _L_ _2_ _9_ 

2 W y 9X X ’ 3X2 ’ 8X3 ’ 0X4 ’3X5’ 


9X 6 ’3X 7 ’3X 8 ’3X 9 


)J2DEHR 


(C2) 


where the partial derivative operators are defined as 



Equation (Cl) is a general expression for any star transit time where the matrix E is 
the only matrix dependent upon 9 . The problem arising from a zero or small 0 may 
be illustrated by examination of this matrix E. The matrix E may be written as 
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E(<M,*) = + *) + sin 0B(( fi,xp) + (1 - cos 9)C {<f>,\p) 


where 


A = 


cos(i^ + <p) 
-sin(\p + <p) 
0 


sin(^/ + <p ) 
cos(ip + <p) 

0 


B = 


C = 


0 0 sin \p 

0 0 cos \p 

sin < p -cos (f> 0 

sin cp sin ip -cos <f> sin ip 0 

sin cp cos xp -cos <p cos xp 0 

-1 


0 0 
From equation (C4) it is clear that for zero 9, 

E(<p,0,xp) = A (<p + xp) 

and for 9 small, for which sin 9 and (1 - cos 9) are small, 


(C4) 


8E ~ 9E 
9 ^0 ~ 0 

and 

3E S 8E Kt _9E_ 3E 

dip 9<j> 9 ^o d ^0 

The matrix multiplier of AX in equation (Cl) therefore has two pairs of columns which 
are nearly equal. Since the matrix [F 2 (x<\] T [F 2 (x«>) k ] <* equation (26) Is determined 

by use of this matrix multiplier, the [F 2 (x«>) k ] T [F 2 (x«>) k ] matrix is clearly singular for 

zero 9 or nearly singular for very small 9. The physical reason for this problem is 
that for small 0, the angles <p Q and > 0 are poorly defined individually and only their 

sum is well defined. Similarly, only the sum of the rates <p and ip are well define 
for small 9. 

The approach to the singularity problem arising from small cone angles was to find 
a new set of parameters Y rather than X, solve a modified least- squares problem for 
AY, and then perform a transformation from the AY vector to the AX vector. During 
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the course of writing a digital computer program for the least- squares equations, 
Mr. Charles B. Grosch of Control Data Corporation suggested the new AY vector: 



(C5) 


which was successfully applied to the singularity problem. Substituting equation (C5) 
into equation (Cl) results in 


0 « J 2 DEHR + VAY 
where the matrix is defined as 

J 2 D(cos 0B + sin 0C)HR 

J 2 D^A^ + sin 6 B^ + (1 - cos 0)C^HR 

tJ 2 D^A^ + sin 0 Bip + (1 - cos 0)C^HR 

J2DEH^,R 


VT 


J2DEHqR 

J 2 d(cos - B , j ,) + sin |(C 4> - C^HR 
1 W 2 d(cos - B„,) + sin Zfo - 

\ J2D € iEHR 

\ J 2 D e2 EHR 



(C6) 


(C 7) 


The subscripts of equation (C7) denote partial differentiation with respect to the particular 
subscripted parameter. Each star transit time then produces an equation of the form 
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given by equation (C6) and the system of resulting equations does not possess a singular- 
ity for small 0. After a solution is made for AY by using the least-squares approach 
described previously, the correction vector AX is found by 



' Ay 2 - Ay 6 
2 sin 0/2 
Ay 3 - Ay ? 
2 sin 0/2 

A y6 

2 sin 0/2 
Ay 7 

2 sin 0/2 


(C8) 


The basis for convergence in the least-squares solution is still the condition that 
AX^ = M as discussed earlier. 
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TABLE I.- THE SEQUENCE OF ROTATIONS AND THE TRANSFORMATION 
MATRICES FOR THE TEN BASIC REFERENCE FRAMES 
USED TO DETERMINE SPACECRAFT ATTITUDE 




TABLE I.- THE SEQUENCE OF ROTATIONS AND THE TRANSFORMATION 
MATRICES FOR THE TEN BASIC REFERENCE FRAMES USED 
TO DETERMINE SPACECRAFT ATTITUDE - Concluded 
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TABLE H.- A PRIORI DATA 



Flight 1 

Flight 2 

Preflight measurement 

y v > deg 

2.8659 

2.8612 

rs» deg 

-0.1967 

-0.2027 

Py, deg 

0 

0 

%, deg 

43.1513 

43.1160 

Ij, slug-ft 2 

130 

135 

l3> slug-ft 2 

8.3 

8.8 

Launch data and flight estimates 

^est* deg/sec 

306.3 

239.6 

ei* deg 

Unknown 

Unknown 

£2> deg 

Unknown 

Unknown 

0 , deg 

75.42 

137.5 

O, deg 

50.95 

58.5 

*0, deg/sec 

19.5 

15.52 

0 O > deg 

Unknown 

Unknown 

0, deg/sec 

286.8 

224.08 

deg 

Unknown 

Unknown 

0, deg 

Unknown 

Unknown 


*0 was estimated with the assump- 


tion that cos 0*1. 




TABLE m.- INDIVIDUAL ELEMENTS OF JDEH MATRIX 


J = 


j COS y 
-sin y cos p 
\sin y sin p 


sin y 

cos y cos p 
-cos y sin p 


0 \ 
sin p 
cos p / 



jcos 62 

sin 

sin 

e 2 

-cos ej sin e 2 \ 

D = 

0 

cos 

e l 


sin ej 


^sin 62 

-sin 

cos 

e 2 

cos ej cos 62 J 


/cos if/ cos <p - cos 0 sin <p sin if/ cos if/ sin cp + cos 9 cos <p sin if/ 

-sin \p cos <p - cos 9 sin cp cos i p -sin if/ sin <p + cos 6 cos cp cos if/ 

\ sin 6 sin cp _ s i n q cos ^ 


H = 


f cos $ 
-cos O sin $ 
sin O sin $ 


sin $ 

cos G cos $ 
-sin 0 cos 4> 


0 

sin 0 
cos 0 


sin if/ sin 9\ 
cos if/ sin 6 
cos 9 j 
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TABLE V.- IDENTIFIED STARS OF PROJECT SCANNER FLIGHT 1 


Boss catalog 
number 

Name 

Bayer name 

Apparent 

visual 

magnitude 

Spectral 

class 

*3584 

Acamar 

6 Eridani 

3.06 

A2 

6274 

Cursa 

/3 Eridani 

2.92 

A3 

6668 

Bellatrix 

y Orionis 

1.70 

B2 

6847 

Mintaka 

6 Orionis 

2.48 

09 

8208 

Tejat Posterior 

p Geminorum 

3.19 

M3 

12407 

Talitha 

l Ursae Majoris 

3.12 

A4 

15145 

Merak 

/3 Ursae Majoris 

2.44 

A0 

16268 

Phekda 

y Ursae Majoris 

2.54 

A0 

17518 

Alioth 

e Ursae Majoris 

1.68 

A0 

*18133 

Mizar 

£ Ursae Majoris 

2.17 

A2 

18643 

Alkaid 

77 Ursae Majoris 

1.91 

B3 

20947 

Alphecca 

a Corona Borealis 

2.31 

A0 

22193 

Kornephoros 

/3 Herculis 

2.81 

G 8 

25180 

Kaus Borealis 

X Sagittarii 

2.94 

K0 

25661 


cp Sagittarii 

3.30 

B 8 

25941 

Nunki 

a Sagittarii 

2.14 

B3 

*26161 

Ascella 

£ Sagittarii 

2.71 

A2 

30942 

Alnair 

a Gruis 

2.16 

B5 


^Double star. 


TABLE VI.- STARS USED TO DETERMINE ATTITUDE OF SCANNER VEHICLES 


Boss catalog 
number 

Name 

Bayer name 

Apparent 

visual 

magnitude 

Spectral 

class 

Scanner flight 1 

*3584 

Acamar 

9 Eridani 

3.06 

A2 

6274 

Cursa 

/3 Eridani 

2.92 

A3 

6668 

Bellatrix 

y Orionis 

1.70 

B2 

6847 

Mintaka 

e Ursae Majoris 

1.68 

A0 

15145 

Merak 

j3 Ursae Majoris 

2.44 

A0 

16268 

Phekda 

y Ursae Majoris 

2.54 

A0 

17518 

Alioth 

e Ursae Majoris 

1.68 

A0 

18643 

Alkaid 

77 Ursae Majoris 

1.91 

B3 

20947 

Alphecca 

a Coronae Borealis 

2.31 

A0 

22193 

Kornephoros 

/3 Herculis 

2.81 

G 8 

*26161 

Ascella 

£ Sagittarii 

_ . 

2.71 

A2 

Scanner flight 2 

519 

Ankaa 

a Phoenicis 

2.44 

G5 

9188 

Adara 

e Canis Majoris 

1.63 

B1 

9443 

Wezen 

6 Canis Majoris 

1.98 

G3 

9886 

Aludra 

77 Canis Majoris 

2.43 

B5 

13926 

Regulus 

Leonis 

1.34 

B7 

*14177 

Algieba 

y Leonis 

2.61 

K 0 

18643 

Alkaid 

77 Ursae Majoris 

1.91 

B3 

25466 

Vega 

a Lyrae 

0.14 

A1 

30491 

Deneb Algiedi 

6 Capricorni 

2.98 

A5 

32000 

Fomalhaut 

a Piscis Austrini 

1.29 

A2 


*Double star. 
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Figure 3 
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Illustration of reticle and pulse groups generated by stars 1 and 2. Lines in 


reticle sketch represent transparent slits. 
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Figure 4.- Star-mapper data-processing procedure. 












Train of coded star presence and absence pulses. 



Cross correlation of star-mapper pulse train with reference code group. 


Star-transit-time detection occurs at 
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Vehicle spin axis 



Figure 10.- Simplified scan of the star-mapper extended field of view about the celestial sphere for the case of zero cone angle. 
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Launch time, vehicle orientation 
and reticle configuration 



Figure 12 .- Flow diagram of error-time-history generation. 





deviation of error in optical- axis -pointing direction, ^(tJ^jdeg 



Figure 13.- 
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Increasing spin rate 
Assumed constant spin rate 



Figure 14.- increase of total angular spin rate during the second data period of Scanner flight 1. 
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Figure 17.- Pointing direction of kg axis for one precession cycle based on the set of nine vehicle motion parameters selected from the 
second data period of Project Scanner flight 1 for time beginning 6h 25m 46.7s UT on August 16, 1966. 








